Issue1575 add pv model #1766
Issue1575 add pv model #1766
Conversation
|
@mwetter All change requests were adressed. Please re-review. |
The proper way to load files is with Modelica.Utilities.Files.loadResource
There was a problem hiding this comment.
@LauraMaier : Pleae see the inline comment for the required changes. I still need to run the model (which fail in Dymola 2024x and Optimica, see inline comment for the cause).
| The Lambert W function solves mathematical equations in which the unknown is both inside and outside of an exponential function or a logarithm. | ||
| </p> | ||
| <p> | ||
| This function is a simple approximation for Lambert W function following Baetzelis, 2016: |
There was a problem hiding this comment.
The citation for Baetzelis is missing.
| info= | ||
| "<html> | ||
| <p> | ||
| The model computes the air mass, which is the number of particles in the atmosphere.</br> |
There was a problem hiding this comment.
The air mass is not equal to the number of particles in the atmosphere.
| "<html> | ||
| <p> | ||
| The model computes the air mass, which is the number of particles in the atmosphere.</br> | ||
| It is based on an exact empirical approach by Kasten et al. and bases on the zenith angle of the object as well as its height. |
There was a problem hiding this comment.
The approach is either empirical (which is here the case) or exact, but not both.
IBPSA/Electrical/BaseClasses/PV/BaseClasses/PVOptical/AirMassModifier.mo
Outdated
Show resolved
Hide resolved
| airMasMod =if (b0 + b1*(airMas^1) + b2*(airMas^2) + b3*(airMas^3) + b4*(airMas^4)) <= | ||
| 0 then 0 else b0 + b1*(airMas^1) + b2*(airMas^2) + b3*(airMas^3) + b4*( |
There was a problem hiding this comment.
Remove the exponent ^1, and use regularization which avoids a time or state event, and it makes the model having a continuous derivative (which is important for NLP).
| HGloHor=realPassThroughHGloHor.y; | ||
|
|
||
| solHouAng = realPassThroughSolHouAng.y; | ||
| solDec=realPassThroughSolDec.y; | ||
| cloTim=realPassThroughCloTim.y; | ||
|
|
||
| k_t =if HGloHor <= 0.01 | ||
| then 0 | ||
| else min(1, max(0, (HGloHor/(GSC*(1 + 0.033* | ||
| cos(360*(Modelica.Constants.pi/180)*cloTim/24/60/60/365)* | ||
| (cos(lat)*cos(solDec)*cos(solHouAng) + sin(lat)*sin(solDec))))))) | ||
| "Factor needed for Erbs diffuse fraction relation"; | ||
|
|
||
| // Erbs diffuse fraction relation | ||
| HDifHor = if HGloHor <= 0.01 | ||
| then 0 | ||
| elseif k_t <= 0.22 | ||
| then (HGloHor)*(1.0 - 0.09*k_t) | ||
| elseif k_t > 0.8 | ||
| then (HGloHor)*0.165 | ||
| else (HGloHor)* | ||
| (0.9511 - 0.1604*k_t + 4.388*k_t^2 - 16.638*k_t^3 + 12.336*k_t^4); |
There was a problem hiding this comment.
Make these assignments in the variable declaration so that we don't mix connect statements with assignments, and use regularization (also as it is only an example, regularization is less critical here).
There was a problem hiding this comment.
Also upload the svg file so that changes can be made.
| parameter Integer nPar "Number of parallel cell circuits on the PV panel" | ||
| annotation(Dialog(group="General")); | ||
| parameter Modelica.Units.SI.Area ACel | ||
| "Area of a single cell. If not found in data sheet, use ACel = ((V_mp0*I_mp0)/(1000*eta_0))/n_ser" |
There was a problem hiding this comment.
Comment is not clear, where should V_mp0 (and other parameters) be obtained? Also, I think you mean nSer.
Maybe such instructions would be better in the info section.
| extends IBPSA.Electrical.Data.PV.Generic; | ||
|
|
||
| parameter Modelica.Units.SI.Efficiency eta0 | ||
| "Efficiency under standard conditions. If not found in data sheet, use eta_0 = ((V_mp0*I_mp0)/(1000*A_cel*n_ser))" |
There was a problem hiding this comment.
Same as in above file, and use the right variable names.
| parameter Modelica.Units.SI.Voltage VMP0 | ||
| "MPP voltage under standard conditions" | ||
| annotation (Dialog(group="Cell specific: Electrical characteristics")); | ||
| parameter Real C |
There was a problem hiding this comment.
Shouldn't this have a unit?
This closes #1575
We added a new PV model to the IBPSA library.
The model has the following characteristics: